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Abstract 

The paper obtains the general form of the cross-covariance function of vector frac- 
tional Brownian motion with correlated components having different self-similarity in- 
dices. 
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1 Introduction 

A p— variate stochastic process X = {X(t) = (Xi(t), ■ ■ ■ , X p (t)),t G R} is said operator 
self-similar (os-s) if there exists a p x p matrix H (called the exponent of X) such that for 
any A > 0, 

X(\t) = fdd X H X(t), (1.1) 

where =fdd means equality of finite-dimensional distributions, and the p x p matrix A^ 
is defined by the power series A^ = e^ logA = H k (log X) k /kl. Os-s processes were 

studied in Laha and Rohatgi (1981), Hudson and Mason (1982), Maejima and Hudson 
(1984), Sato (1991) and other papers. A Gaussian os-s process with stationary increments 
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is called operator fractional Brownian motion (ofBm). For p = 1 the class of ofBm coincides 
with fundamental class of fractional Brownian motions (fBm) (see e.g. Samorodnitsky and 
Taqqu (1994)). Recall that a fBm with exponent H € (0, 1) can be alternatively defined 
as a stochastically continuous Gaussian process X = {X(t),t € M} with zero mean and 
covariance 

2 

EX(s)X(t) = y (|s| 2 ^ + \t\ 2H -\t- s\ 2H ), t,s£R, (1.2) 

where a 2 = EA 2 (1). The form of covariance of general ofBm seems to be unknown and 
may be quite complicated. The structure of ofBm and stochastic integral representations 
were studied in Didier and Pipiras (2008). 

A particular case of os-s processes corresponds to diagonal matrix H = diag(i?i, • • • , H„). 
In this case, relation (jl.ip becomes 

(X 1 (\t),---,X p (\t)) = fdd (A Hl Xi(t),---,A^X p (t)). (1.3) 

Below, a p— variate process X satisfying (|1.3j) for any A > will be called vector self-similar 
(vs-s) and a stochastically continuous Gaussian vs-s process with stationary increments (si) 
will be called a vector fractional Brownian motion (vfBm). 

Vs-s processes (in particularly, vfBm) seem to be most useful for applications and sta- 
tistical analysis of multiple time series. They arise as limits of normalized partial sums of 
multivariate long memory processes with discrete time, see Marinucci and Robinson (2000), 
Davidson and de Jong (2000), Chung (2002), Davidson and Hashimadze (2008), Robinson 
(2008), Lavancier et al. (2009). In particularly, vs-s processes appear in the OLS estima- 
tion in multiple linear regression model (Chung, 2002), the multiple local Whittle estimation 
(Robinson, 2008) and the two-sample testing for comparison of long memory parameters 
(Lavancier et al., 2009). 

Note from (|1.3|) that each component X{ = {Xi(t),t € M}, i = 1, • • • ,p of vs-s process is 
a (scalar) self-similar process, the fact which is not true for general os-s processes. 

The present paper obtains the general form of the (cross-)covariance function of vs-s si 
process X with finite variance and exponent H = diag(iii, • • • , H p ), < Hi < 1. According 
to Theorem 12.11 below, under some regularity condition, for any i, j = 1, • • • ,p, i ^ j with 
Hi + Hj ^ 1, there exist c^-, Cji E K such that for any s, t £ M 



°»°J f„. r,,M„ _!_ ,...,/, i, i //,+//, , ; _ _ ,,//, + // 

where a 2 := var(Aj(l)) and 



cov(A,( S ),A,(t)) = -Y{ Cij ( s )\s\ H ^+c ji (t)\t\ H ^ - Cji (t-s)\t-s\ H *+ H i}, (1.4) 



f ca, t > 0, 

c « (t) : = U. t <0. M 

A similar expression (involving additional logarithmic terms) for the covariance cov(Aj(s), 
Xj(t)) is obtained in the case Hi + Hj = 1. We prove Theorem 12.11 by deriving from vs- 
s si property a functional equation for the cross-covariance of the type studied in Aczel 
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(1966) and Aczel and Hosszu (1965) and using the result of the last paper (Theorem 12.21 
below) about the uniqueness of this equation. Section [3] discusses the existence of vfBm 
with covariance as in (jl.4p . We start with a double sided stochastic integral representation 
similar to Didier and Pipiras (2008): 

X(t)= [ {((t-x)^-- 5 -(-x)^-- 5 )A + + ((t-x) H ~- 5 -(-x) H -- 5 )A\w(dx), 

(1.6) 

where H — .5 := diag(-ffi — .5, ■■■,H P — .5),x+ := max(x,0),x_ := max(— x, 0), A+, A- 
are real p x p matrices and W{dx) = (Wi(dz), • • • , W p (dx)) is a Gaussian white noise with 
zero mean, independent components and covariance EWi(dx)Wj(dx) = Sijdx. According 
to Proposition 3.1, if < Hi < 1, + H 3 i ^ 1, i,j = 1, ■ ■ • ,p then the cross-covariance of 
X in (jl.6p has the form as in (jl.4p with 

= 2cij4>ij/(Ji(jj, (pij := B(Hi + .5,Hj + .5)/sm((Hi + Hj)ir), (1.7) 

where the matrix C = (c^ ) is given by 

C := cos(H-k)A + A* + + A^A*_ cos(-Htt) (1.8) 
- sm(Hir)A + A*_ cos(Hir) - cos(Hir)A + A*_ sm(Hir). 

Here and below, A* denotes the transposed matrix, sin(.£f7r) := diag(sin(//i7r), • • • , sm(H p ir)), 
cos(Hir) := diag(cos(-fTi7r), • • • , cos(// p 7r)). 

2 The form of the covariance function of vs-s process 

Recall that a random process X = {X(t),t G M} has stationary increments (si) if {X{t + T) 
-X(T),t£R} = fdd {X(t) - X(0),t G M} for any T G M. 

Theorem 2.1 Let X = {X(t),t G M} 6e a 2nd order process with values in W. Assume 
that X has stationary increments, zero mean, X(0) = 0, and that X is vector self-similar 
with exponent H = diag(iii, • • • , H p ), < H j < 1 (i = 1, • • • ,p). 

Moreover, assume also that for any i,j = l,---,p, the function t i— > ~EXi(t)Xj(l) is 
continuously differentiable on (0,1) U (l,oo). Let a\ > denote the variance of Xi(l), i = 
1, • • • 

(i) If i =j, then for any (s,t) G M 2 , we have 

2 

EXi{s)Xi(t) = ^- {\s\ 2m + \t\ 2 ^ - \t - s\ 2m } . (2.1) 

(ii) If i ytz j and Hi + Hj / 1, then there exists Cij,cji G R such that for any (s,t) G M 2 , 
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(Hi) If i jtz j and Hi + Hj = 1, then there exists dij,fij G R such that for any (s,t) G M 2 , 
we have 

EXi(s)Xj(t) = {t%(|s| + |t| - |s - t|) + /^(ilog |t| - slog |s| - (t - s) log |t - a|)} . 

(2.2) 

(w;J TTie matrix R = (Rij)i,j=i,—,p is positive definite, where 

{1, i=j, 
Cjj + Cji, i ^ j, -ffj + Hj ^ 1, 
dij, i / j, Hi + Hj = 1. 

Proof, (i) Follows from the well-known characterization of covariance of (scalar-valued) 
self-similar stationary increment process. 

(ii)-(iii) From part (i), it follows that EpQ(t) - X;(s)) 2 = of |i - s\ 2H > and hence t £ 

R} is stochastically continuous on the real line, for any i = 1, ■ ■ ■ ,p. Whence, it follows 
that FiXi(s)Xj(t) is jointly continuous on R 2 and vanishes for t = or s = 0, for any 

hj = - ,p- 

Denote r(s,i) := EXj^X^i), # := ± (#* + #,•). Let R+ := {u : u > 0}, R_ := {u : u < 
0}. From vs-s and si properties we obtain 

r(Xs,Xt) = X 2H r(s,t), (2.3) 
r(s,t) = r{s + T,t + T)-r(s + T,T)-r{T,t + T)+r{T,T) (2.4) 

for any reals s, t, T and any A > 0. Substituting s = t = u, A = l/\u\ into (|2.3h one obtains 

r(u,u) = 2k±|u[ 2H , uGR±, 2k± := r(±l,±l). (2.5) 

Substituting s = t = -l,T = 1 into flU]) and using r(0,0) = r(l,0) = r(0, 1) = yields 
r(— 1, —1) = r(l, 1). Next, substituting T = u,s = t = v — u into (|2.4p and using (|2.5p with 
k := K + = k_ one obtains 

r(u, u) + r(w, u) = r(u, u) + r(u, u) — r(v — u, v — u) 

= 2k(\v\ 2H + \u\ 2H -\u- v\ 2H ) (u, v G R). (2.6) 

Next, let 

g±(t) := r(±l,t) (i G R). 

Then (|2.3p implies 

r(s,t) = |s| 2H c/ ± (±t/s) (s G R±, t G R). (2.7) 
Equation (j!T4"|) with s = 1 and (|27f|) yield, for all t G R 



./-;-(" - r/ -nr'"., + (^ 



)-^ + (^) 



and for all i £ R, equation (|2.4p with s 
g-(t) = (T-1) 2H 

- (T-1) 2H 



#4 



.94 



-1 and (1221) yield 

T + 1 \ _ - r 2ll 

T - 1 
T 



T - 1 

We claim that the general solution of functional equations ([2? 
• When H / 1/2, 



+ T^V(1) (T>1). (2.9) 
2.91) has the form: 



9+(t) 



' c' + c'\t\ 2H -c'\l-t\ 2H , 
c' + c\t\ 2H -c'\l-t\ 2H , 



c 1 - t 



2# 



MO 



L C + c\t\ 2H 
c + c'\t\ 2H -c'\t + l\ 2H , 
c + c'\t\ 2H -c\t + l\ 2H , 
[c + c\t\ 2H -c\t + l\ 2H , 



t < 0, 

< t < 1, , 

t > 1, 

t < -1, 
-1< t < 0, 
f > 0; 



(2.10) 



(2.11) 



• When H =1/2, 
9-(t) 



r/(tiog|t|-(t-i)iog|t-i|), t<o, 

\d(lAt) + /(tlog|t|-(i-l)log|t-l|), t>0, 

f d(l A |t|) + /(tlog|t| - (t + l)log|i+l|) , t < 0, 

\/(ilog|i|-(i + l)log|i + l|), t>0, 



(2.12) 
(2.13) 



with some c,c',d,f G R. 

It follows from ([XT0>([27l3l and ([277)1 that the covariance r(s,t) 



(s,i) G R 2 has the form as in ([Q|> . (f272"|) . with Cy = 2c/ciiO- 



= EXi(s)X,-(t) for 
, — 2c' /ui(Jj in the case 
2// = Hi + Hj y^z l i and = 2d/aicrj, /y = 2/ /<Ji&j in the case 2// = Hi + Hj = 1. 

To show the above claim, note, by direct verification, that (|2. 10|) - (|2.13l) solve equations 
(p^]) - p^|) . Therefore it suffices to show that p7lQ|l - f[2Tl3| ) is a unique solution of (|278j) - (|23jh 
Let t > 1. Differentiating (|2.8p with respect to £ leads to 



g' + (t) = (T+l) 



Let x := t, y 
(1, oo) x (0, oo) 



(T + t)/T. Then (x,y) G (1, oo) 2 and the mapping (£, T) 
(1, oo) 2 is a bijection. Equation (|2.14p can be rewritten as 



x + y — 1 



2H-1 



where 



F(x,y) := 



<?4» + 



x + y — 1 



2H-1 



g'+{y). 



(2.14) 

(x,y) : 

(2.15) 
(2.16) 



x + y — 1 

Equation (|2.15p belongs to the class of functional equations treated in Aczel (1966) and 
Aczel and Hosszu (1965). For reader's convenience, we present the result from Aczel (1966) 
which will be used below. 
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Theorem 2.2 (Aczel and Hosszu (1965)) There exists at most one continuous function f 
satisfying the functional equation 

f(F(x,y)) = L(f(x),f(y),x,y) (2.17) 

for all x, y £ (A, B) and the initial conditions 

f( ai ) = bi, f(a 2 ) = b 2 (ai,a 2 € (A, B), a x ^ a 2 ) 

i/-F is continuous in (A,B) x (A,B), L(u,v,x,y) is strictly monotonic in u or v and F is 
intern (i.e., F(x,y) E (x,y) for all x ^ y € (A,B)), where (A,B) is a closed, half-closed, 
open, finite or infinite interval. 



Unfortunately, equation (|2,15p does not satisfy the conditions of Theorem 12,21 (with / = 
<7+ and (A, B) = (1, oo)), since F in (12. 16ft is not intern. Following Aczel and Hosszu (1965), 
we first apply some transformations of (|2.15p so that Theorem 12.21 can be used. 

Note, taking T = t/(t — 1) in ([231]) one obtains 

g' + (t) = K(t)g' + (£—^, (2.18) 



where 

(2t - l) 2H - 1 

t 2H-l + {t _ 1 



K ® ■= ^V/Zi^-i - ( 2 - 19 ) 



Let F(x) := F(x,x) = 2 x-i ■ Then F is strictly increasing from (l,oo) onto (l,oo). Let 

G(x,y) := F {F(x,y)) = ■ , 2.20 

x + y — 1 



with F 1 (y) = y + y/y(y - 1). Note F(x,y) > 1, G(x,y) > 1 for (x,y) £ (l,oo) 2 . From 
(12.181) with t = G(x, y) one obtains 

g' + (G(x,y)) = g' + (F(x,y))K(G(x,y)). (2.21) 

Combining (12. 15H and (12.21H one gets 

g' + (G(x,y)) = L(g' + (x),g' + (y),x,y), (2.22) 

where 

L(u,v,x,y) := -) -) v ) K(G(x,y)). (2.23) 

+ y — 1/ \x + y — 1/ / 

The fact that G in (|2.20p is intern follows from its definition and monotonicity in x and y, 
implying 

x = G(x,x) < G(x,y) < G(y,y) = y 



for any x < y, x,y G (l,oo). Since L is monotonic in u or v, so Theorem 12,21 applies to 
functional equation (|2,22p and therefore this equation has a unique continuous solution g' + 
on the interval (1, oo), given boundary conditions g' + (ai) = hi (i = 1, 2), 1 < a 2 < ai < oo. 
Form of g + on [1, +oo) : 

• Assume H ^ 1/2. Let ai := 2,b x := cH(2 2H - 2). In view of (f2T8l) . the other 
boundary condition can be defined by a 2 = a?/(2ai-l) = 4/3, 6 2 := g' + {2)/K{2) and 
so equations (12.22D and (|2.14p have a unique solution for single boundary condition 
g' + (2) = cH(2 2H - 2). (See also Aczel and Hosszu (1965, p.51).) Since g' + (t) = 
2Hc(t 2H ~ 1 — (t — l) 2H ~ l ) is a solution of (|2.14j) with this boundary condition, it 
follows that this solution is unique. Hence it also follows that 

g + (t) = c' + ct 2H - c{t - l) 2H ie[l,oo), (2.24) 

for some c' € M. 

• When H = 1/2. A particular solution of (|2. 14[) when t > 1 is log(i) — log(i — 1). For 
the same reason as above, the general solution of (|2.14p is thus d'(log(t) — log(t — 1)) 
where d! £ BL It follows that 

g + (t) = d + d! {t\ogt - (t - l)\og{t - I)) tG[l,oo), (2.25) 

for some d, d' G M. 

Form of <?+ on (0, 1) : 

Putting t = 1 in (|2.8p results in 

(T+1) 2 V(^)+T 2 V(^) = <7 + (l)[(T+ir + r 2 *-l]. 
Whence, for s := T/(T + 1) S (0, 1) and using (|2.24p one obtains 

9+ (s) = - S 2H g + (l/ S ) + g + (l)[l + s 2H -(l-s) 2H ]. 

(2.26) 

This gives, for s € (0, 1), 

• when H / 1/2, g + (s) = c' + cs 2H - c'(l - s) 2H , 

• when H = 1/2, g+(s) = ds + d'(slogs — (s — l)log(s — 1)). 
Therefore, relations (I2,10p and ()2.12p have been proved when t > 0. 

Form of g+ on (— oo, 0) : This case follows from (|2.8|) taking T = —t. 
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Form of g_ : The relations ([2"TT]) and ([2TT3J) are deduced from (|23U|) and ([27T21) thanks 
to relation (j2l)j) . 

(iv) Follows from the fact that R is the covariance matrix of random vector {Xx(l)/cr\, 
■■■,X p (l)/a p ). 

Theorem 12.11 is proved. 



3 Stochastic integral representation of vfBm 

In this section we derive the covariance function of vfBm X = {X(t),t £ M} given by 
double-sided stochastic integral representation in (11.61) . Denote 

v P P p 

a tj + ■= Yl a ik a i> a if : = S a ib a jk> a tf : = 2 °^ a 7fc' % + := Yl a ik a tk> 

k=l k=l k=l k=l 

where A + = [afA , A_ = ( a~- ) are the p x p matrices in f)l .6j) . Clearly, 



A + A\ = [at 3 + ), A-At = [ar-y A+A*_ = (a±~ ) , A^A* + =(a- j 

Note, each of the processes Xi = {Xi(t),t £ M} in (jl.6p is a well-defined fractional Brownian 
motion with index Hi £ (0, 1); see e.g. Samorodnitsky and Taqqu (1994). 

Proposition 3.1 The covariance of the process defined in hl.b}) satisfies the following prop- 
erties 

(i) For any i = 1, • • • ,p the variance of Xi{l) is 

2 B{Hj + .h, Hj + .5) f ++ . + , 

^ = sin(gfTr) ^ +a " -^(^K^ /• 

fw,) //-ffj + Hj ^ 1 i/ien /or any s,f£t, t/ie cross- covariance ~EXi(s)Xj(t) of the process 
in hi. 6)) is given by ( |i.^| ), wii/i 

o-iO-j __ ^(ffj + .5,Hj + .5) 
2 Qj : ~ sm((Hi + Hj)7r) 

- cos(if i7 r) + ay" cos(i^-Tr) - a±" sin((i7; + #,)tt)} . (3.2) 



x < aj 3 



(Hi) If Hi + Hj = 1 then for any s,f£l, i/ie cross- covariance EXi(s)Xj(t) of the process 
in hi. 0(1 is given by h2.2\) . with 



in 

OiOjdij := B(Hi + .5, Hj + .5) 

( sm(HiTr) + sm(Hjir) , , , , ,1 . „. 

x| 2 L ^(«5 + + %")-«5"-% + )> (3-3) 

ao-jhj := (H j -H i )(a++ -arr). (3.4) 
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Remark 3.2 Let Hi + Hj 7^ 1, i, j = 1, • • • ,p and let Qj, 0jj be denned as in (jl.7p . 

^From Proposition 3.1 (|3.ip . (|3,2p it follows that the matrix C = (cjj) satisfies f)1.8|) . In 
this context, a natural question arises to find easily verifiable conditions on the matrices C 
and H such that there exist matrices A + , A- satisfying the quadratic matrix equation in 
(jl.8p . In other words, for which C and H there exists a vfBm X with cross-covariance as 
in (fTIj) ? 

While the last question does not seem easy, it becomes much simpler if we restrict the 
class of vfBm's X in (jl.6p to causal representations with = 0. In this case, equation 
(|1.8p becomes 

(5 = cos(^)^+^;. 

Clearly, the last factorization is possible if and only if the matrix cos{H)~ l C is symmetric 
and positive definite. 



Proof. For all s, let 



I+(s) := I (( s -x)^- 5 -(-x)^- 5 )^(dx), 

' fl *- 5 - (-x)5 1- - 5 ) Wfe(dx). 



s — X 



Using the above notation, Jf;(s) = £fe=i (a^i^(a) + a ik I ik (s)) and 

EX,( S )X,(i) = a± + E/+( S )// 1 (t) + a+-E/+( S )/r 1 (t) 



(3.5) 



Let i^j + Hj 7^ 1. From Stoev and Taqqu (2006, Th. 4.1), taking there a + = 1, a = 0, 
H(s) = Hi and H(t) = Hj, we obtain 



EI+(s)I+(t) = 4> H 



cos (( Hj -H { f-- (g ' + 2 ^ )7r sign( S )) | a |*+*« 

cos - + {Hi+ 2 Hj)ir *m(tj) \t\ H ^ 



- cos {{Hj-Hi) 



7T (JTj + fl--)7T 



signfs -t) lis - i|- ff »+- H J 



where 

IpH 

Therefore, 

e4 (*) 



r(fli + .5)r(flj + .5)r(2 -H { - Hj) B(Hi + .5, Hj + .5) 



n(Hi + Hj)(l - Hi - Hj 
B(Hi + .5,Hj + .5) 



sm((Hi + i^-) 71 ") 



sin((ifi + Hj)-K) 



bij(s)\s\ H > +H i + 6 i< (*)|t| Hi+H ' - Ma " " 
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bij(s) 



where 

cos(iIj7r), if s > 
cos(Hjir), if s < 0. 

Similarly, taking a + = 0, a - = 1, H(s) = Hi and H(t) = Hj in Stoev and Taqqu (2006, 
Th. 4.1), we obtain 



Finally, 
E^(*Kji(f) 



g(gj + .5,gj + .5) 
mn((fli + Hj-)tt) 



Hi+Hj 



l {s>t} 



\s - x) Hi - 5 (x - t) H '- 5 dx - 1 {S>0} f\s - x) Hi - 5 x H i-- 5 dx 



-l {t<0} j\-x)^-\x-t) H ^dx) 
B(H { + .5, Hj + .5) ((a - t) J i+H ^ - s+ i+Hi - (3.6) 



and 



EJ^( S )/+(t) = 5(fr i + .5,fl> + .5)f(t-a 



\Hi+Hj Hi+Hj Hi+Hj 



Substituting these formulas into (|3.5p we obtain (|3.2p and (|3.ip . 

Next, let Hi + Hj = 1. We get similarly from Theorem 4.1 of Stoev and Taqqu (2006) 

1 



(3.7) 



7T 



B{Hi + .5,Hj + .5) 



^sm(Hi7r)(\s\ + \t\-\s-t\) 



■ cos (-£/j7r) (slog |s| — tlog |t| — (s — i) log |s — t\) 



and 



7T 



+ .5, Hj + .5) 



7T 



sin(^7r)(|s| + |t| - \s - t\) 



+ cos (-ffj7r) (slog |s| — t log |i| — (s — t) log |s — t\) 

Expressions (|3.6p and (|3.7p remain true when i% + Hj = 1 and they can be rewritten as 

1 



E/+(s)/- 1 (i) = EI a ( s )I+(t) 



+ .5, #j + .5)(|s| + \t\ - |s - t|). 



Therefore, using (|3.5p . we obtain (|3.3p and (|3.4p . Proposition 3.1 is proved. 
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